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Abstract. Testing of various classes of life distributions has been ad- 
dressed in the literature for more than forty-five years. In this paper, we 
consider the problem of testing exponentiality (which essentially implies 
no aging) against positive aging which is captured by the fairly large 
class of New Better than Used in Expectation (NBUE) distributions. 
These tests of exponentiality against NBUE alternatives are discussed 
and compared. The empirical size of the tests is obtained by simula- 
tion. Power comparisons for different popular alternatives are done us- 
ing Monte Carlo simulation. These comparisons are made both for small 
and large sample sizes. The paper concludes with a discussion in which 
suggestions are made regarding the choices of the test when a particular 
alternative is suspected. 
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1 Introduction 

The exponential distribution is characterized by the lack of memory prop- 
erty; a constant hazard rate and a constant mean residual life function. 
This is the most commonly used distribution in reliability and survival 
analysis, primarily because of its mathematical simplicity. The assump- 
tion of exponential life times essentially means that a used item is stochas- 
tically as good as a new one; hence there is no need to replace a unit which 
is working. However, this is not always a realistic assumption; and age 
does have an effect on the residual life time. Positive (negative) aging 
means that age has an adverse (beneficial) effect, in some probabilistic 
sense, on the residual life. Hence, it is of interest to check possible depar- 
ture from exponentiality in the data. 

A hierarchy of aging classes has been proposed in the literature to 
model the effect of aging. One such aging class is the so-called New Bet- 
ter than used in Expectation (NBUE) class, which is defined below. 



Definition 1 ; A non-negative random variable X with finite mean ^ is 
said to possess the NBUE property ifepit) = Ef{X — t\X > t) ^ Ep{X) 
i.e 



Important results for this class are discussed in Barlow and Proschan 
(1975). Specifically the NBUE family of life distributions is important in 
the study of replacement policies. In fact, Marshall and Proschan (1972) 
have proved that the average waiting time between any two consecutive 
failures when no planned replacement policies is adopted is smaller than 
or equal to the similar quantity when an age replacement policy is adopted 
if and only if the life distribution is NBUE. 

Furthermore, this average waiting time is the same under these two 
policies only when the system life is exponentially distributed. Hence, if 
the average waiting time between the consecutive failure is an important 
criterion in deciding whether to adopt an age replacement policy over the 
failure replacement policy for a given system, then a reasonable way to 
decide would be to test whether the life distribution of the given system 
is exponential. Rejection of the exponential hypothesis on the basis of 
the given data would imply as favoring the adoption of age replacement 
policy. 

Doksum and Yandel (1984) and Spurrier (1984) were the first one to 
review the tests of exponentiality. This was followed up with the works of 
Ascher (1990) and Henze and Meintanis (2005). Anis and Dutta (2010) 
discuss tests of exponentiality against IFR alternatives; while Anis (2010) 
reviews the procedures of testing exponentiality against DMRL alterna- 
tives. We contribute to the literature by examining the tests of exponen- 
tiality against NBUE alternatives. 

This paper is organized as follows. We introduce the different test 
procedures in Section 2. The results of a Monte Carlo study on the size 
of the tests is reported in Section 3. We present in Section 4 the power 
of these tests against the commonly used alternatives to the exponential 
model. Finally concluding remarks are made in Section 5. 




2 The Test Statistics 



The problem of testing exponentiality against NBUE alternatives is more 
than three decades old. Different types of approaches are used in deriving 
the test statistics. The first test procedure was proposed by Hollander 



and Proschan (1975) based on a distance measure, de Souza Borges et al. 
(1984) proposed a test based on the coefficient of variation. Belzunce et 
al. (2000) had considered the same problem by means of the dispersion of 
residual lives; while Belzunce et al. (2001) proposed a test based on the 
right spread order. Fernandez-Ponce et al. (1996) suggested a test based 
on the right-spread function to characterize different partial orderings 
between life time distributions. Anis and Mitra (2011) have generalized 
the Hollander-Proschan statistic to obtain a family of test statistics. 
Formally we are concerned with testing the following hypothesis: 

Ho : F G v.s Hi : F G NBUE 

based on a random sample Xi , X2 , . . . , of size n from an absolutely 
continuous distribution function F{x) = 1 — e~^^,x ^ and A > 0, 
typically unknown. Let F = 1 — F he the corresponding survival function. 
Let = X(o) ^ X(^i^ ^ X(2) ^ . . . ^ -^{n) the corresponding order 
statistics. Let X denote the sample mean. 
We shall now briefly discuss these tests. 

2.1 The First Test Procedure 

Hollander and Proschan (1975) suggested a test statistic by considering 
the integral 

poo 

7(F) = / F{t){eFiO)-eFit))dFit) (1) 
Jo 

as a measure of deviation from Hq towards the alternative Hi . The sample 
counterpart to 7(F) is obtained by substituting the empirical distribution 
function F„ for F in (1) above. On simplification, they obtained 

i=i ^ ■' 
where X(^j^ denotes the j-th order statistic. In order to make the test 
statistic scale invariant, they suggest using Ti = K/X. They have proved 
that their test is consistent and asymptotically normal. The null hypoth- 
esis is rejected for large value of Ti. 

2.2 Test due to Koul (1978) 

By definition it follows that F is NBUE if 
ry 

/ F{x)dx ^ lJiF{y) < y < 00 



where // = F{y)dy < oo. Then define 

D:= sup \- rF{x)dx- F{y)\ . (3) 

0<y<oo I A* Jo J 

Observe that if F is exponential, then D = 0. The farther D is from 
0, the more evidence there is for F E Hi. Substituting for F by the 
empirical distribution function Fn and simplifying Koul (1978) obtained 
the statistic 

T2 = Dn = D{Fn)= max [wni--\ (4) 



where 

™ E-=i(n-J + l){X(,)-X(,„,)}- 

The test rejects Hq in favour of //i if T2 is large and is consistent. The 
null distribution of this test procedure is derived in Barlow and Doksum 
(1972). 



2.3 Test based on the Coefficient of Variation 

de Souza Borges et al. (1984) introduced a simple test based on the sample 
coefficient of variation. They considered the parameter 

A{F)= / F{t){tx-eF{t)}dt (5) 
Jo 

as a measure of the deviation from the null hypothesis of exponentiality. 
They proved the following characterization of exponentiality for life dis- 
tributions which are continuous NBUE. 



Theorem : Let F he a continuous life distribution and NBUE. A neces- 
sary and sufficient condition for F to he exponential is that CV{F) = 1, 
where CV{F) denotes the coefficient of variation of F 

Using the above characterization theorem, they suggest a test based on 
the sample coefficient of variation S/ X, where 

n 

n ^-^ ' 



Under iJo, SjX is asymptotically normal with mean 1 and variance 1/n. 
A large sample level a test rejects in favour of the NBUE property 



if Ta = ^/n{S/X — 1) < — Zq,, the standard normal a-quantile. (The test 
rejects Hq in favour of the NWUE property if ^/n{S/X — 1) > Za)- They 
also showed that the test is consistent. 



2.4 Test due to Aly (1990) 

Aly considered the parameter 



7F = / (1 - y) {1 + log(l - y)} dF~\y) ^ y ^ 1 
{ei.(0) - eF{t)} dF^t) 



as a measure of NBUE-ness. Substituting for the empirical distribution 
function and simplifying he obtained 

. . E (^) } (^) {M^} 

(6) 

He has proved that y/nT^^ is asymptotically standard normal; and is 
consistent. The test rejects the nul hypothesis of exponentiality at (ap- 
proximate) level a if 

y/^{^{Fn)-\nXn] ^ 

— ^ Zi^a 



<7nXn 



where 



An=i+-x:iog fi-^v 

i=i ^ ^ 



and zi-Q, is the (1 — a) quantile of the standard normal distribution. 



2.5 Test based on a Quantile Dispersion Measure 

Fernandez-Ponce et al. (1996) proposed the test statistic T5 based on the 
Right Spread (RS) function. The RS function of the random variable X 
is defined as 

poo 

S+{u) = E[{X-F^\u))+)= Fx{x)dx (7) 

Jf-\u) 



where x+ = max{x, 0} and F ^(x) is the left continuous inverse (quantile 
function), Qx{u) = F~'^{x) = mf{x : Fx{x) ^ u}; Vn G (0, 1). 

Definition 2 : F is more new better than used in expectation than G, 
NBUE 

(denoted by F < G) if 

NBUE 

It can be shown than F is NBUE if and only if F < G , where 

Their test is motivated by considering the following integral as a mea- 
sure of deviation, for a given F^, from Hq to Hi : 



Jo 



(8) 



where L^+(0) is the length of the segment determined by the secant to 
S^{u) from U = to U = u in the x-axis. Making the change u = Fx{t) 
one has 

nF) = I {l-Ljf (0)}rfF(t). (9) 

F(t) 

The parameter ^{F) can be viewed as follows. Define D{t) = 1 — (0). 
Now, D{t) = Vt G M 44> Hq is true. To show this equivalence observe 
that, 

a^(F) = i?^(l-LjW(0)). (10) 

So '^{F) is simply an average value of the deviation from Hq towards Hi. 
Note that ^{F) ^0 V F G U and ^{F) = if and only if F G Hq. 

The sample analogue of the parameter 'F{F) is the basis for their 
proposed test statistic. Given a random sample of size n, they obtain 
the test statistic T5, by replacing F{t) by its natural estimates, namely 
the empirical distribution function. Thus Fernandez-Ponce et al. (1996) 
obtain 

( i\2l ,(n-k + l)(X(.)-X(. l^] 
Tr,=^(F„) = l--y l-^^i^ '] itzlVl (11) 



with the assumption that X^gj = 0. 



Fernandez-Ponce et al. (1996) state that is asymptotically unbiased 
and converges almost surely to ^{F). Though the exact distribution of 
could not be found, the critical points for sample size n = 3 (1) 20 (5) 35 
for different significance levels have been tabulated in Fernandez-Ponce 
et al. (1996). 

2.6 Test based on Dispersion of Residual Life 

Belzunce et al. (2000) proposed a new test for Hq based on dispersion 
of residual life. Before describing their test, we need to give some back- 
ground definitions. 

Definition 3 -.(Hickey (1986)) A random variable Y is more dispersed 
in dilation than a random variable X (denoted by X Y) if 

E {<p{X - E{X))) ^ E {<P{Y - E{Y))) (12) 

for every convex function (j). 

Belzunce et al. (1997) have proved that 

X is NBUE if and only if Xt >du ^, Vt > (13) 

where Xt = {X — t\X > t} is the remaining life at age t and Y is 
an exponential random variable with mean equal to E{X), i.e. E{Y) = 
E{X). 

Using the characterization in (13) above, Belzunce et al. (2000) con- 
sider 

^NBUE {X) = [ Flit) {AduiXt, Y)} dFxis) (14) 
as a measure of departure from Hq to Hi, where 

roo r-oo 

Ada (X, Y) = / FUt)dt - Fix) - / (t)dt + E{Y) (15) 

Jlx Jly 

with Ix and ly being the smallest values of the supports of F^ and Fy, 
which are assumed to be finite; and where Y is an exponential distribution 
with mean equal to E{X). 

Now, Ambue{X) = under the null hypothesis of exponentiality, 
and under Hi, Ai^bue{X) > 0. Based on a random sample of size n, the 
empirical analogue is given by 

Anbue {n)^-Y^ n{n - i)^ A„(i) + (16) 

i=0 ^ ^ 



where 



A„(i) 




1 



(17) 



n 



(18) 



a=i+l 



n — 2a + i + 1. 



(19) 



The nuh hypothesis of exponentiahty is rejected for large values of 
^NBUE{n). Since ANBUE{n) is not scale invariant, Belzunce et al. (2000) 
suggest using 



They have derived the exact distribution of Tg and have shown that Tg 
is asymptotically normal. Under Hq, the limiting distribution of V AbnTg 
is standard normal. They have proved that the test is consistent and 
the critical points based on the exact distribution is obtained for n = 
2 (1) 20 (5) 60. For large values of n, the null hypothesis of exponentiahty 
is rejected if V 45nTg > Zq,. 

2.7 Test based on the Right Spread Order 

Belzunce et al. (2001) suggested a new test procedure based on the right 
spread order. 

Recall that the Right Spread Function has already been defined in 
sub-section 2.5 above. It should be noted that the right spread function 
can be considered as a measure of spread to the right of every quantile 



A random variable X is said to be less than Y in the right spread 
order (denoted by X Y) if 



E{{X- F-i(n))+} ^E{{Y- FyHu))+} Vn G (0, 1). (21) 



It is interesting to note that the right spread function of a random 
variable X is also related to its mean residual life function epit) = E{X — 
t\X > t) by the relationship 



Te = l\*^BUE[n) 



^NBUEin) 



(20) 



Qx{u). 



E{{X- E^\u))+} = (1 - u)eHF^Hu)). 



(22) 



Then it is easy to show that 



X e NBUE <^ X i^Rs Y 



(23) 



where Y is an exponentially distributed random variable such that E{Y) = 
E{X). 



Belzunce et al. (2001) use the above characterization in (23) above to 
construct a family of tests for NBUE aging class. In fact their proposed 
test statistic is similar to their test statistic A^^(X, y) proposed to test 
the right spread order; i.e. 

H'^ ■.X=RsY v.s. H[ : X ^rs Y. 

More specifically for testing the null hypothesis of exponentiality against 
NBUE alternatives, they suggest taking A^^(X, y) where Y is an expo- 
nential random variable with E{Y) = E{X), as a measure of departure 
from Hq towards Hi. Then it can be shown that 



1 poo 

/\%BjjE^l^ls{X,Y) = -E{X){l-a){2-a)- Up) Fx{x)dx dp 

o Jo Jf^\p) 

(24) 

where 

fp(i-l) ^^p^a 
a ^ p ^1. 



Ja{p) 




Then under the null hypothesis of exponentiality l\^^u^{X) = 0, but 
under ffi, A'^^^^{X) > 0. Based on a random sample, the sample ana- 
logue of the measure in (25) above is given by 



A^W^i,--,^n) = ^X(l-a)(2-a)-^f:{L«(i)-J„(l)(l ' ^ 



and 



n 
(25) 



where ^ a, ^ > a. 



Since this statistic is not scale invariant, they suggest using 



Belzunce et al. (2001) have found the exact distribution of Tj, under 
the nuh hypothesis of exponentiahty. Further, they have proved that T^ 
is asymptoticaUy normal and under Hq, the Hmiting distribution of 

Thus the null hypothesis of exponentiahty is rejected at level a if for 
large n, -^^^j'j^^^^Tj > Za- They have also proved that the test is 
consistent. 



2.8 Test based on Moment Inequality derived from comparing 
Life with its Equilibrium Form 

It should be noted that associated with the life distribution X is the no- 
tion of "random remaining life" at age t, denoted by Xf. It is easy to 
show that Xf has a survival function Ft{x) = F{x + t)/ F{t), x,t ^ 0. It is 
well known that, see for example Ross (2003), that Xt converges weakly 
to a non- negative random variable X with survival function Wpix) = 
j^CF{u)du,x^O. 

In an interesting paper Mugdadi and Ahmad (2005) showed that the 

~ St 

notion that X is NBUE is equivalent to the notion that X ^ X; where 

~ St 

X denotes the equilibrium life, and ^ denotes less than or equal to in the 
usual stochastic order. 

Then by Theorem 3 of Mugdadi and Ahmad (2005), Hq can be tested 

~ ~ St 

by comparing X and the equilibrium life X. Since X ^ X means that, 

/"OO 

/ F{u)du ^ uF{x) (27) 

J X 

^ (^j^F{u)du^dF{x)^^ (28) 

^ rF\x)dx ^ ^ (29) 
Jo 2 

Hence as a measure of departure from Hq, we may consider 

^8 = - !^^^-E{mmiX,,X2))] . (30) 



Then Tg can be estimated by 

ft4(;i^EE{f-M^<.^]^ (3.) 

Using U-statistics theory, Mugdadi and Ahmad (2005) show that ^/n(Tg — 
Tg) is asymptoticahy normal with mean and variance 

^2 = Var + I - 2 £ xdF{x) - 2XiF(Xi)^ . (32) 

Then under Ho,a^ = 1/12. The test rejects Hq for large value of the 
statistic Tg. 

2.9 A Generalization of the HoUander-Proschan Type Test 

Anis and Mitra (2011) proposed a family of test statistics by generalizing 
the Hollander and Proschan (1975) approach. Essentially they consider 
the measure 

roo 

7,(F) = / {my {epiO) - epit)} dF{t) (33) 
Jo 

as a measure of departure from Hq towards the alternative Hi . The sample 
analogue of ^j{F) reduces to 

Clearly ^j{Fn) is an L-statistic, being a linear combination of order statis- 
tics. As 7j(F„) is not scale invariant, Anis and Mitra (2011) suggest us- 
ing To = ^*j{Fn) = ^j(Fn)lyi,, where = E{X). The null hypothesis is 
rejected for large values of 7*(F„). They proved that 7*(F„) is asymp- 
totically normal and the test is consistent. Based on the sample observa- 
tion, the scale invariant statistic 7*(F„) is estimated by Tq = 7*(F„) = 

7;(F„)/X;, where X„ = E7=l^^/n■ 

Very recently, Anis and Basu (2011) have obtained the exact sampling 
distribution of 7j'(F„). They have also given a set of critical values for the 
commonly used significance levels and sample sizes n = 2 (1) 25 (5) 100. 
They conclude from their table that the convergence to normality is ex- 
tremely slow. 



3 Size of Different Tests 



In a large simulation study, we calculate the empirical size of the tests 
under consideration; i.e, the rejection probability under Hq. We do this for 
different sample sizes and for nominal size (level of significance) of 5%. To 
calculate the empirical size of the tests we simulate n observations from 
an exponential distribution and compute the respective test statistic and 
check whether this particular realization of the test statistic accepts or 
rejects the null hypothesis of exponentiality. This procedure is repeated 
for 10^ times and the proportion of times that the test statistic rejects 
the null hypothesis is observed. Thus we estimate the size for the different 
tests. The standard errors of the estimated sizes are less than or equal to 
y^(l25/T¥ = 0.0016. Ah simulations were done using MATLAB 7.8 on 
PC platform. 

3.1 Size for small sample 

Hollander and Proschan (1975) refers to Barlow (1972) for the small 
sample critical values. Fernandez-Ponce et al. (1996) and more recently 
Belzunce et al. (2000) provide the critical values for small sample sizes. 
Anis and Basu (2011) have provided the critical values for the Anis and 
Mitra (2011) test Tq. Anis and Basu (2011) have also provided the critical 
values for the Hollander and Proschan (1975) test statistic based on the 
exact distribution which they have obtained. Table 1 shows the simulated 
sizes for these four tests for n = 5 (1) 15 and for a nominal significance 
level a = 5%. Also included are the simulated sizes for the Anis-Mitra 
test To for j = 1. It is seen that the Anis-Mitra test with j = 0.25 appears 
to be very conservative being always less than 0.05. The test T5 due to 
Fernandez-Ponce et al. (1996) overshoots the nominal size. 

3.2 Size for moderate sample size 

We again compare the above mentioned four tests procedures for mod- 
erate sample sizes. We consider n = 16 (1) 20 (5) 30. Similar trends are 
seen. The Anis-Mitra (2011) test Tq is most conservative; while T5 due 
to Fernandez-Ponce et al. (1996) overshoots the nominal size most of the 
time, as seen in Table 2. 

3.3 Size for large sample size 

For large sample sizes, most of the test statistics are normally distributed. 
However, whenever the exact critical values are available, we have used the 



exact value. The test proposed by Fernandez-Ponce et al. (1996) has a very 
comphcated statistic and its asymptotic distribution has not been derived. 
Hence this test procedure has not been considered for large sample. 

From Table 3, we see that the test T2 due to Koul (1978) and the test 
T3 based on the coefficient of variation due to de Souza Borges et al. (1984) 
are very conservative. The test T7 based on the right spread order given by 
Belzunce et al. (2001) is also conservative. The test Tg based on dispersion 
of residual lives constructed by Belzunce et al. (2000) overshoots the 
nominal size for most of the sample sizes considered. The test T4 due to 
Aly (1990) performs poorly and overshoots the nominal size for all sample 
sizes. The test Tg based on the equilibrium distribution given by Mugdadi 
and Ahmad (2005) also exceeds the nominal size. Finally the tests Tq 
given by Anis-Mitra (2011) and Hollander Proschan (1975) Ti perform in 
a mixed manner; for some sample sizes the tests are conservative while 
sometime they overshoot the nominal size. 

4 Power of the Tests 

The power of the tests under comparison has been evaluated by Monte 
Carlo simulation. The power estimates were calculated as the proportion 
of 10^ Monte Carlo samples that resulted in rejection of the null hypoth- 
esis of exponentiality at significance level a of 5% for the alternative dis- 
tributions considered. The standard errors of the estimated probabilities 
are bounded above by y^O. 25/10^ = 0.0016. As mentioned in the previous 
section, the simulations were done using MATLAB 7.8 on PC platform. 
The power of the tests is simulated from the following distributions: 

— The Weibull distribution with density 6x^~^exp{x^),9 ^ l,x ^ 0, 
denoted by W{9); 

— The Gamma distribution with density {-r(0)}~^x^~^ea;p(— x), ^ l,x^ 
0, denoted by r(6'); 

— The linear failure rate distribution with density (1 + 9x)exp{—x — 
9x^/2), e ^ 0, X ^ 0, denoted by LFR{9). 

Note that these distributions are generally used as alternatives to the 
exponential model. It is easy to show that the Weibull and gamma distri- 
butions belong to the NBUE class for 9 > 1; while the LFR distribution 
is NBUE for ^ > 0. The null distribtuion is obtained at = 1 for the first 
two; whereas the LFR distribution reduces to the exponential distribution 
at 9o = 0. 

In a large simulation study, we simulated the power for different al- 
ternatives for sample sizes n = 5 (5) 100. However, for brevity we shall 



report our findings for n = 5 (5) 30 (10) 100. We simulated the pow- 
ers for Weibull alternatives for shape parameter 9 = 1.10 (0.10) 1.50. 
For the gamma distribution, the shape parameter 9 selected were 9 = 
1.20 (0.20) 2.00. Finally, for the LFR distribution, the shape parameter 
9 chosen were 9 = 0.25 (0.25) 1.25. It should be noted that for these pa- 
rameter values the distributions are NBUE. Hence a "good" test should 
enjoy high power. 

It may be pointed out that Anis and Mitra (2011) have observed 
an error in the computational formula for the Hollander-Proschan test 
statistic. Anis and Basu (2011) have made a comparative study of the 
Hollander-Proschan statistic Ti and the Anis-Mitra statistic Tq with j = 
1. It is further reinforced from the simulation results in Table 4-6, that 
we should use the correct statistic as given by Anis and Mitra (2011). 
Though the difference is of the order of 1 /n? , it is appreciable for small 
sample sizes, as expected. 

4.1 Power for small sample size 

We are able to compare only four test procedure when the sample size is 
small. These tests are due to Hollander and Proschan(1975) Ti; Fernandez- 
Ponce et al. (1996) T5; Belzunce et al. (2000) Tq and Anis and Mitra 
(2011) Tq. For small sample comparative study we chose the sample size 
n = 5 (5) 25. Tables 4-6 show the simulated power of these tests against 
the Weibull, Gamma and the LFR alternatives. We observe from these 
tables that the performance of each of these tests is poor when the NBUE 
property is marginal. This is expected because the test can easily confuse 
and consider the observations to have come from an exponential distri- 
bution. 

From Tables 4 and 5 it is observed that if the Weibull or Gamma 
alternatives are suspected then it best to use the Anis-Mitra test Tq for 
very small {n ^ 10) sample size; while it is preferable to use the test T5 due 
to Fernandez-Ponce et al. (1996) for moderate sample size (15 ^ n ^ 25). 
It is also observed that it is better to use the Anis-Mitra test with j = 1 
instead of the original Hollander-Proschan test Ti. The test Tq due to 
Belzunce et al. (2000) comes a poor third. 

For the LFR alternatives, it is observed from Table 6, that the Anis- 
Mitra test Tq with j = 0.25 performs best for all sample sizes n ^ 25. The 
performance of test T^ due to Fernandez-Ponce et al. (1996) and Tq due to 
Belzunce et al. (2000) based on dispersion of residual lives is comparable. 



4.2 Power for large sample size 



We have compared the performance of eight of the tests against WeibuU, 
Gamma and LFR alternatives for different values of the shape parameter 
6 and for sample size n = 30 (5) 100. For brevity of space we give in Table 
7-9 the power comparisons for sample sizes n = 30, 40, 50, 75 and 100 for 
different alternatives. The test due to Fernandez-Ponce et al. (1996) 
has not been considered as neither the distribution of the test statistic is 
derived nor the critical values (based on simulation) are provided. 

It is seen from Table 7 and 8 that for Weibull and gamma alternatives, 
the Anis-Mitra test Tq with j = 1 performs best. The performance of 
the original Hollander-Proschan test Ti is comparable, though the power 
figures are smaller (in view of the inherent error in the test statistic) . The 
performance of the Anis-Mitra test Tq with j = 0.25 is also very good. 
For Weibull alternative, the test due to Mugdadi and Ahmad (2005) Tg 
is worst. For gamma alternatives, the worst performance is exhibited by 
the tests due to Aly (1990) and Mugdadi and Ahmad (2005) Tg. As an 
example, when the alternative is Weibull with = 1.5, with a moderately 
large enough sample of size n = 30, the Anis-Mitra test with j = 1 is 
able to take the correct decision 81.87% of the time. The corresponding 
figure for the Mugdadi and Ahmad test Tg is mere 30.73% only. Even 
with a large sample of size n = 100, the Mugdadi and Ahmad test Tg 
is able to take the correct decision only 77.28% of the time when the 
scale parameter is 9 = 1.3, the corresponding figures for the Hollander- 
Proschan test Ti is 93.17%, while for the Anis-Mitra test Tq with j = 1 
is 95.44% and for j = 0.25, it is 92.84%. 

Similarly when the Gamma parameter ^ = 2, (a moderate NBUE 
property) the test T4 due to Aly (1990) is able to take correct decision 
only 1% of the time with a sample size of n = 30. In contrast the test 
Tg due to Mugdadi and Ahmad (2005) has a success rate of 24.39%. For 
these parameter values, the test Tq with j = 1 shows the estimated power 
as 79.04. 

For the LFR alternatives, it is seen that the best performance is ex- 
hibited by the test T4 due to Aly(1990). It is a clear winner. The worst 
performance is shown by Tg due to Mugdadi and Ahmad (2005). 

General Remarks : It is seen that the test Tq-i due to Anis and Mitra 
(2011) with j = 1 out performs the test Ti due to Hollander and Proschan 
for all sample sizes and for all alternatives considered. 



5 Concluding Remarks 



In this work the performance of several tests for exponentiahty against 
NBUE alternatives has been studied for different sample sizes. We shall 
summarize the major findings. 

It is better to use the test To{j = 1) due to Anis and Mitra (2011) 
instead of the asymptotically equivalent test Ti due to Hollander and 
Proschan (1975). 

For small and moderate sample sizes, the test ro(j = 0.25), To(j = 1) 
and Ti are very conservative. For large sample sizes the tests T2 and 
are very conservative. On the other hand test T^^Tq and Tg overshoots 
the minimal size. 

For very small sample size (n ^ 10) it is better to use ro(j = 0.25) due 
to Anis and Mitra (2011) if the alternative is suspected to be Weibull or 
Gamma. If the sample size is small (15 ^ n ^ 25), then the performance 
of the test due to Fernandez-Ponce et al. (1996) is best and is recom- 
mended for use. However, if the sample size is small and the suspected 
alternative is LFR, then the best performance is exhibited by the test 
To(j = 0.25) due to Anis and Mitra (2011). 

If the sample size is large (n ^ 30) and the alternative is suspected to 
be Weibull or gamma alternatives, then it is recommended that TQ{j = 1) 
due to Anis and Mitra (2011) be used. The worst performance in this case 
is exhibited by T4 due to Aly (1990) and Tg due to Mugdadi and Ahmad 
(2005). In contrast if the alternative is suspected to be LFR, then the best 
test to use is T4 due to Aly (1990). Even in this case the performance of 
the test Tg due to Mugdadi and Ahmad (2005) is worst. Hence the test 
proposed by Mugdadi and Ahmad (2005) should be avoided. 
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Table 4. Estimated power for the Weibull alternative for small sample size 





oaiiipic 


Anis and Mitra (2011) 


llOiidliLiei alUA 


J7 (3i lidliLie^ JT uii(_,e 




vveiULiiii (7 ) 


Size 








r^iosciian iiy(O) 


ei ai. I ly yu ) 


ei ai. I zuuu f 




Tl 


3 — 


J — U.OU 


i — 1 


1 


-Lb 


-te 







7 on 
( .uu 


qq 


7 

( .uo 


/ .11 


7Q 


7 Oft 




1 n 

lU 


8.51 


8.63 


8.70 


ft ftO 
O.OZ 


O.oD 


ft 

O.oD 


9 = 1.10 


15 


10.01 


10.15 


10.29 


8.88 


12.2 


9.54 




20 


11.12 


11.21 


11.35 


10.56 


12.64 


11.22 




25 


11.73 


12.02 


12.29 


11.71 


17.26 


11.34 




5 


9.39 


9.48 


9.42 


9.59 


8.14 


8.43 




10 


12.68 


12.74 


12.83 


13.09 


9.71 


12.67 


6 = 1.20 


15 


16.50 


16.50 


17.07 


16.47 


20.56 


15.59 




20 


21.05 


21.45 


22.34 


19.96 


24.18 


18.95 




25 


23.06 


23.71 


24.65 


23.74 


33.29 


21.86 




5 


11.73 


11.95 


11.98 


12.19 


11.38 


11.83 




10 


18.53 


18.54 


18.9 


19.14 


14.76 


17.68 


6 = 1.30 


15 


26.92 


27.66 


28.38 


25.42 


32.43 


23.64 




20 


31.90 


33.01 


34.68 


33.15 


37.90 


29.59 




25 


38.21 


39.66 


41.36 


40.16 


51.58 


36.20 




5 


13.84 


13.88 


14.13 


14.32 


13.92 


14.00 




10 


25.90 


26.24 


27.28 


27.67 


21.03 


24.25 


e = 1.40 


15 


35.88 


36.82 


38.05 


37.58 


44.14 


33.59 




20 


46.62 


48.06 


50.29 


48.52 


53.46 


44.16 




25 


55.63 


57.97 


60.00 


58.01 


67.87 


52.29 




5 


17.52 


17.86 


18.04 


18.28 


16.98 


16.72 




10 


33.61 


34.44 


35.72 


36.14 


28.87 


31.79 


e = 1.50 


15 


47.59 


49.50 


51.06 


48.60 


57.02 


44.48 




20 


59.77 


61.60 


63.94 


62.78 


68.89 


57.98 




25 


70.49 


72.23 


74.13 


73.17 


81.84 


66.87 



Table 5. Estimated power for the Gamma alternative for small sample size 
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Table 6. Estimated power for the LFR alternative for small sample size 
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Table 7. Estimated power for the WeibuU alternative for large sample size 
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81.87 


81.58 


58.39 


68.59 


83.49 


75.35 


77.39 


30.73 




40 


89.37 


90.84 


92.06 


91.91 


75.21 


83.33 


91.89 


87.82 


89.42 


57.60 


e = 1.5 


50 


94.57 


95.63 


96.53 


96.24 


84.96 


90.79 


96.35 


93.40 


94.34 


75.69 




75 


99.57 


99.70 


99.83 


99.64 


97.00 


98.51 


99.55 


99.46 


99.22 


95.58 




100 


99.93 


99.95 


99.98 


99.96 


99.63 


99.77 


99.93 


99.91 


99.89 


99.42 



Table 8. Estimated power for the Gamma alternative for large sample size 







Anis and Mitra 


Hollander &: 


Koul 


do Souza 


Aly 


Belzunce 


Belzunce 


Mudgadi 


Gamma 


Sample 




(2011) 




Proschan 


(1978) 


Borges et 


(1990) 


et al. 


et al. 


& Ahmed 




Size 




ToU) 




(1975) 




al. (1984) 




(2000) 


(2001) 


(2005) 




n 


j = 0.25 


j = 0.50 


i = i 


Ti 


Ta 


Ts 


Ti 




T7 


Ts 




30 


14.50 


14.89 


15.67 


14.39 


7.46 


8.64 


6.38 


12.86 


13.32 


1.32 




40 


15.72 


16.51 


17.33 


16.52 


9.28 


10.43 


7.26 


14.30 


14.98 


2.49 


e = 1.2 


50 


18.14 


19.13 


20.72 


20.04 


11.67 


12.97 


8.77 


17.22 


18.62 


4.38 




75 


27.41 


29.63 


33.10 


25.95 


16.37 


16.80 


11.12 


26.22 


23.33 


8.31 




100 


32.16 


34.26 


38.13 


31.52 


20.99 


22.04 


14.75 


30.37 


28.39 


13.41 




30 


26.95 


28.57 


31.11 


29.75 


16.32 


18.45 


5.37 


25.07 


26.68 


3.64 




40 


33.07 


35.22 


38.48 


38.31 


22.67 


25.02 


6.95 


31.11 


34.09 


8.33 


e = 1.4 


50 


37.26 


40.24 


44.52 


44.40 


28.26 


30.11 


9.40 


36.19 


39.93 


13.94 




75 


56.62 


60.78 


66.78 


58.25 


42.01 


41.98 


15.91 


54.20 


52.98 


26.59 




100 


68.56 


72.91 


78.40 


72.59 


55.17 


53.96 


23.36 


65.79 


66.63 


39.61 




30 


42.27 


44.77 


49.04 


48.29 


28.20 


32.08 


3.95 


39.71 


42.71 


8.45 




40 


52.05 


55.48 


60.83 


59.98 


39.33 


42.62 


6.06 


49.50 


54.25 


18.30 


e = 1.6 


50 


59.84 


63.92 


69.77 


68.61 


49.36 


51.06 


8.60 


57.54 


62.90 


28.84 




75 


82.09 


86.10 


90.19 


86.15 


70.99 


68.87 


18.92 


79.69 


81.16 


52.36 




100 


89.64 


92.81 


95.79 


93.43 


84.39 


79.66 


32.11 


88.17 


89.87 


68.77 




30 


57.28 


60.53 


65.60 


64.49 


40.72 


46.42 


2.14 


54.21 


58.70 


15.72 




40 


68.78 


72.96 


78.04 


78.45 


57.74 


60.30 


4.17 


67.11 


72.24 


32.58 


e = 1.8 


50 


78.61 


82.66 


87.13 


86.35 


70.81 


69.79 


6.68 


75.24 


81.36 


48.11 




75 


93.59 


95.98 


97.88 


96.32 


89.35 


85.14 


19.70 


92.13 


93.66 


74.16 




100 


97.72 


98.76 


99.54 


99.11 


96.49 


93.17 


35.02 


97.11 


98.17 


87.50 




30 


70.68 


74.26 


79.04 


78.67 


55.26 


59.35 


1.07 


66.47 


71.81 


24.39 




40 


81.85 


85.50 


89.63 


89.51 


72.96 


74.35 


2.31 


79.48 


84.63 


47.24 


e = 2.0 


50 


89.37 


92.27 


95.07 


94.91 


84.29 


83.04 


4.99 


87.00 


91.81 


64.06 




75 


98.16 


99.00 


99.70 


99.19 


96.65 


94.27 


17.31 


97.76 


98.25 


87.98 




100 


99.70 


99.87 


99.96 


99.95 


99.37 


98.27 


36.28 


99.64 


99.80 


96.47 



Table 9. Estimated power for the LFR alternative for large sample size 







Anis 


and Mitra 


Hollander & 


Koul 


dc Souza 


Aly 


Belzunce 


Belzunce 


Mudgadi 


Lit H 


Sample 




(2011) 




Proschan 


(1978) 


Borges et 


(1990) 


et al. 


et al. 


& Ahmed 


{9) 


Size 




ToU) 




(1975) 




al. (1984) 




(2000) 


(2001) 


(2005) 




n 


j = 0.25 


j = 0.50 


i = i 


Ti 


Ta 


73 


T4 


Te 


Tt 


Tg 




30 


14.52 


14.68 


14.49 


14.03 


7.12 


9.20 


27.05 


15.23 


13.74 


1.42 




40 


18.28 


17.99 


17.75 


17.37 


9.77 


12.45 


30.96 


18.44 


17.03 


3.01 


e = 0.25 


50 


22.02 


21.68 


20.82 


20.45 


12.60 


16.09 


35.55 


21.59 


19.96 


5.40 




75 


35.99 


36.11 


35.06 


27.30 


18.27 


23.86 


44.89 


35.07 


26.97 


12.33 




100 


42.70 


42.43 


40.70 


33.39 


23.84 


31.78 


52.82 


41.59 


32.69 


19.67 




30 


25.14 


24.90 


24.37 


23.30 


12.59 


16.42 


48.94 


24.73 


22.96 


2.65 




40 


30.84 


30.79 


30.29 


31.08 


18.47 


24.54 


58.11 


32.86 


30.57 


7.40 


e = 0.50 


50 


39.09 


38.88 


37.14 


36.13 


23.00 


30.38 


64.26 


38.31 


35.52 


12.60 




75 


60.54 


60.54 


58.65 


50.37 


35.67 


47.29 


77.42 


59.35 


49.07 


29.25 




100 


72.38 


71.60 


69.13 


62.00 


46.72 


60.99 


86.33 


70.47 


60.60 


45.36 




30 


34.24 


33.85 


32.77 


31.22 


17.78 


23.48 


66.20 


33.36 


31.00 


4.70 




40 


44.04 


43.72 


42.32 


40.55 


25.43 


33.35 


74.40 


42.92 


39.97 


11.49 


e = 0.75 


50 


51.21 


51.03 


48.74 


49.47 


32.30 


43.60 


81.43 


51.68 


48.33 


19.77 




75 


75.50 


75.25 


73.53 


65.55 


49.24 


63.37 


91.33 


74.40 


64.14 


44.19 




100 


85.72 


85.41 


83.39 


77.97 


63.36 


77.48 


95.92 


84.82 


76.54 


64.02 




30 


41.50 


41.32 


40.23 


39.06 


22.12 


30.14 


77.66 


40.54 


38.52 


6.75 




40 


52.14 


51.66 


50.06 


49.68 


32.02 


42.43 


85.22 


51.55 


48.53 


16.77 


e = 1.00 


50 


62.41 


62.16 


59.58 


58.91 


41.11 


53.73 


90.12 


61.58 


57.59 


28.36 




75 


84.92 


84.67 


82.81 


76.59 


60.63 


74.31 


96.35 


83.86 


75.48 


56.77 




100 


92.50 


91.97 


90.46 


86.88 


75.08 


86.98 


98.74 


91.64 


85.54 


76.75 




30 


48.14 


48.30 


47.22 


44.58 


27.06 


35.31 


84.82 


46.06 


43.89 


9.01 




40 


59.20 


59.37 


57.48 


56.67 


37.33 


49.08 


91.19 


58.65 


55.36 


20.63 


e = 1.25 


50 


69.38 


69.08 


67.02 


66.81 


48.79 


61.33 


94.23 


69.08 


66.00 


35.62 




75 


90.07 


89.85 


88.40 


83.26 


68.38 


81.94 


98.65 


89.32 


82.32 


66.17 




100 


95.72 


95.41 


94.22 


92.13 


82.82 


92.24 


99.50 


95.09 


91.32 


84.84 



